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o : 

D , Abstract 

o- 

■ In this article we introduce algorithms which compute iterations of Gauss-Manin 
I connections, Picard-Fuchs equations of Abelian integrals and mixed Hodge structure 

■ of afline varieties of dimension n in terms of differential forms. In the case n = I such 
r— I , computations have many applications in differential equations and counting their limit 

• cycles. For n > 3, these computations give us an explicit definition of Hodge cycles. 

<: 

^ ' 1 Introduction 

'■ 

^ I The theory of Abelian integrals which arises in polynomial differential equations of the 
type X = P{x,y), y = Q{x,y) is one of the most fruitful areas which needs a special 

. attention form algebraic geometry. The reader is referred to the articles [3], [H] and 

J> I for a history and applications of such Abelian integrals in differential equations. In this 
article we deal with computational aspects of such integrals. All polynomial objects which 

. we use are defined over C. 

CN I Let us be given a polynomial / in n+1 variables xi, X2, ■ ■ ■ , Xn+i, a polynomial differen- 

' tial n-form u) and a continuous family of n-dimensional oriented cycles 6t C Lt := f~^{t). 

• The protagonist of this article is the integral J^^ oj, called Abelian integral. Computations 
"^H ! related to these integrals becomes easier when we put a certain kind of tameness condition 

' on / (see For such a tame polynomial we can write uj as: 

a: , 



in 



X 
J3 



where 7/^,/3 G / is a class of differential n- forms constructed from a basis of the Milnor 
vector space of / and p^'s are polynomials in t (see ^for the algorithm which produces 
Pps). The Guass-Manin connection Vu has the following basic property 

The above term can be written in the form with pps rational functions in t with 
poles in the critical values of / (see ^JSJfor the algorithm which produces pps). The n-th 
cohomology of a smooth fiber Lt is canonically isomorphic to where is the 

restriction of polynomial differential i-forms to Lt , and carries two natural filtrations called 
the weight and the Hodge filtrations (both together is called the mixed Hodge structure). 

^Math. classification: 14C30, 32S35 
Keywords: Mixed Hodge structures of affine varieties, Gauss-Manin connection, Brieskorn module. 
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These filtrations are generalizations of classical notions of differential forms of the first, 
second and third type for Riemann surfaces in higher dimensional varieties. The reader 
who is not interested in the case n > 1 is invited to follow the article with n = 1 and 
with the usual notions of differential forms of the first, second and third type. How to 
calculate these filtrations by means of differential forms is the main theorem of U and 
related algorithms are explained in fJHl Last but not the least, our protagonist satisfies a 
Picard-Fuchs equation Yli=oPi(^)w ~ ^^^^^ Pi'^ polynomials in t. The algorithm 
which produces pp^s is explained in ^ The theory of Abelian integrals can be studied 
even in the case n = 0, i.e. / is a polynomial in one variable. Since some open problems, 
for instance infinitesimal Hilbert Problem (see 0), can be also stated in this case, we have 
included ^ All the algorithms explained in this article are implemented in the library 
brho.lib of SINGULAR and throughout the article we mention the related procedures of 
brho.lib. Nevertheless, the reader may want to implement the algorithms of this article 
in any other software in commutative algebra. 

The main theorem of ,5j(see also 3Z|) does not give a basis of the Brieskorn module 
compatible with the mixed Hodge structure (see Definition In ^Ifll we obtain such bases 
for some examples of / by modifying the one given in ^ (we do not have a general method 
for every /). Applications of our computations in differential equations and particularly 
in direction of the article 1 is a matter of future work. 

2 Tame polynomials and Brieskorn modules 

We start with a definition. 

Definition 1. A polynomial / € C[x] is called (weighted) tame if there exist natural num- 
bers ai,a2, . . . , On+i S N such that Sing{g) = {0}, where g = /d is the last homogeneous 
piece of / in the graded algebra C[x], deg{xi) = Oi. 

The multiplicative group C* acts on C"'^^ in the following way: 

A* : (xi, X2, . . . , xn+i) ^ (A^^xi, A'^^xa, . . . , A""+ix„+i), A G C* 

The polynomial (resp. the polynomial form) uj in C""^^ is (weighted) homogeneous of 
degree d G N if A*(u;) = A'^w, A G C*. Fix a homo geneous polynomial g of degree d 
and with an isolated singularity at G C^^^. Let Ag be the affine space of all tame 
polynomials / = /o + /i + • • • + fd-i + 9- The space Ag is parameterized by the coefficients 
of fi,i = 0,1, d — 1. The multiplicative group C* acts on Ag by 

A . / = = X-'^fo + A-'^+Vi + • • • + A-Vd + 9 

The action of A G C* takes A • / = biholomorphically to / = 0. 

Let / G Ag. We choose a basis := {x^ | /3 G /} of monomials for the Milnor 
C-vector space 

V := C[x]/iacob{g) 

In Singular one can get x^ using kbase command. Using this command deg(x^) is not 
a decreasing sequence. In brho . lib the procedure okbase makes a permutation on the 
result of kbase and gives us x^'s with deg(x^) decreasing. We will fix the order obtained 



2 



by okbase and use it for many examples throughout the text. Define 

n+l 

(3) Wi:=^, l<i<n + l, r] := {J^i-'^T'^WiXid^i), U := f-\t),t G C 

1=1 

n+l 

:= + l)wi, r]p := x^r], ujp = x'^dx 13 e I 

i=l 

where dxi = dxi A ■ ■ ■ A dxi^i A dxi^i A • • • A dxn+i. Note that Ap = — 1, \i turns 
out that x^ is also a basis of V/ := C[x\ / ]acob{f) and so / and g have the same Milnor 
numbers (see the conclusion after Lemma 4 of jS]). We denote it by /i . We denote by P 
the set of critical points of / and by C := f{P) the set of critical values of /. We will use 
also P for a polynomial in <C.[x\. This will not make any confusion!. 

Let 0*, i = 1, 2, . . . , n + 1 (resp. fi*-, j € N U {0}) be the set of polynomial differential 
z-forms (resp. homogeneous polynomial differential i-forms) in C""*"^. The Milnor vector 
space of g can be rewritten in the form V := j^^^n . The Brieskorn modules 

ttI ttI Tjll ttII 

^ '~ df A + fi!^"-i ' ~ ^ ~ df A dQ'^-^ 

of / are C[t]-modules in a natural way: t.[u)] = [fuj], [uj] G H' resp. € H" . They are 
defined in the case n > 0. The case n = is treated separately in ^ 

Definition 2. Let H be one of H' or H" . If H = H" then by restriction of a; on L^, c G 
C\C we mean the residue of in Lc and by j^uj, 5 G HniLc,"^) we mean Resi{j^). 
It is natural to define the mixed Hodge structure of H as follows: WmH, m € Z (resp. 
F^H, A; € Z) consists of elements uj G H such that the restriction of uj on all Lc, c € C\C 
belongs to W„i7"(Lc,C) (resp. F^H'''{Lc,C)). 

In the case of a tame polynomial we have 

{0} = Wn^i dWnd Wn+i = H, {0} = C F" C • • • C C F° = F 

Each piece of the mixed Hodge structure of is a C[t]-module. In the same way we 
define the mixed Hodge structure of the localization of H over multiplicative subgroups of 
C[t]. In the case n = 1 our definition can be simplified as follows: We have the filtrations 
{0} = WoCWiCW2 = H and = C C F'^ = H , where 

Wi = {lo H \ UJ restricted to a regular fiber has not residue at infinity } 

F^ = {uj £ H \ UJ restricted to a regular fiber has poles of maximum order 1 at infinity} 
In particular we get 

Wi n F^ = {uj H \ u! restricted to a regular comactified fiber is of the first kind} 

For the notion of compactification of and infinity see |T| and . 

The projection of F* in Gr^H := Wm/Wm _i gives us the filtration F* in Gr^H and 
we define Gr'^pGr^H = F^/F^+^. 

Definition 3. Suppose that H is freely generated module. The set B = Um,fceZ-B^ C H is 
a basis of H compatible with the mixed Hodge structure if form a basis ofGrpGr^H. 
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3 Quasi-homogeneous singularities 

Let / = 5 be a weighted homogeneous polynomial with an isolated singularity at origin. We 
explain the algorithm which writes every element of H" of 5 as a C[t] -linear combination 
of LOps. Recall that 

OXj OXi OXi OXj 

where €ij = if i < j and = if i > j and dxi, dxj is dx without dxi and dxj (we have 
not changed the order of dxi, dx2, ... in dx). 

Proposition 1. For a monomial P = x^ we have 

(4) -^-Pdx = — T^gdx + dgA d{\ ^ — ^- ^-XjPdx^,dxj) 

0x4 d.Ag — Ui OXi dAg - ai 

Proof. The proof is a straightforward calculation. 

— -da A dixjPdxi, dxj) 

dAp-at 



-1 y dg d{x,P) ^ dg djxjP) 
dAg — a-i ^ dx-i dxi dxi dx.\ 



1/3 — OXj OXi OXi 

■1 // , dq BP 9(7 

—-adg - - E aM + 1)) 



dAg — a-i dxi dxi dxi 



dAg — ai dxi dxi dxi 

^ 3+1- 



We use the above Proposition to write every Pdx E il"^"^ in the form 



□ 



(5) Pdx = ^pp{g)LOp + dgAd^, p/g G C[t], ^ G n'''\degip(3ig)u;p,dg Ad^) < deg(Pdx) 

/3G/ 

• Input: The homogeneous polynomial g and P G C[x] representing [Pdx] G H" 
Output: P/3,/3 G / and ^ satisfying © 
We write 

(6) Pdx = cpx^.dx + dg Arj, deg{dg A r/) < deg{Pdx) 

Then we apply @ to each monomial component P^ of dg Ar] and then we write 

each §§-dx in the form ©. The degree of the components which make Pdx not to 
be of the form (jSJ always decreases and finally we get the desired form. 
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This algorithm is implemented in the procedure linearl of the library brho .lib. To find 
a similar algorithm for H' we note that if 77 € $7" is written in the form 

(7) r^ = Y,Pp{9)m + d9'\i + di^, pp e C[tli,ii ^ n^-' 

where each piece in the right hand side of the above equality has degree less than deg(?7) 
then 

(8) dr] = ^{pp[g)Ap + p'p{g)g)ujp - dg A 

pel 

and the inverse of the map C[t] — > C[i], p — > Af^.p + p' .t is given by X]i=o^i** ~^ 
X](^=i Af^+i ■ Therefore, we can recover pps using linearl for dr]. The obtained al- 
gorithm is implemented in the procedure linear2 of the library brho. lib. Later we will 
introduce the procedure linear (resp. linearp) which works for arbitrary tame poly- 
nomial and is an extended version of linearl (resp. Iinear2). For this reason, the 
procedures linearl, linear2 are not available to the user. 

Theorem 1. For a weighted homogeneous polynomial g, the set 
with 

Bn+i = {oJp\Ap = n-k + l},B^ = {ujp\n-k<Ap<n-k + l}, 

is a basis of H" compatible with the mixed Hodge structure. The same is true for H' 
replacing top with rjp. 

This theorem is due to J. Steenbrink [7j and its generalization for an arbitrary tame 
polynomial is given in ^ 



4 A basis of H' and H" 

Proposition 2. For every tame polynomial f G Ag the forms ujp, (3^1 (resp. tjp, P I) 
form a basis of the Brieskorn module H" (resp. H' ) of f. More precisely, every uj € $7"^"^ 
(resp. UJ € $7" J can be written 

(9) u; = ^ppif)u;p + df A d^, pp G C[t], ^ S deg(p^) < - Ap 

pel 

(resp. 

(10) u = Y.pp{f)vp + df Ai + d^i, ppe C[t], e G deg(p;3) < - Ap 
) 

This Proposition is proved in f5 Proposition 1. The proof also gives us the following 
algorithm to find all the unknown data in the above equalities. 
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• Input: The tame polynomial / and P E C[x] representing [Pdx] G H" . 
Output: p/3,/3 € / and satisfying © 

We use the algorithm of ^and write an element uj € $7"+^, deg(u;) = m in the form 
u} = ^pp{g)ujp + dg A dtp, pp E C[t], G ^^''"S deg(p/3(5)u;^), deg(d5 Ad^p) <m 

This is possible because g is homogeneous. We have 

^ = + df ^di) + J, J = -P(}{f))^i3 + d{g - f) Adip 

f3£i pal 

The degree of is strictly less than m and so we repeat what we have done at the 
beginning and finally we write cj as a C[t]-linear combination of cj^'s. 

The algorithm for H' is similar. The statement about degrees is the direct consequence 
of the proof and ©. The procedure linear takes a tame polynomial / and P G C[x] 
and returns a list. The first entry is a 1 x /i matrix C = and a x matrix 

representing a (n — l)-form ^. The base ring must contain at least one parameter and C 
is written in the first parameter. The procedure linearp for H' is similar to linear for 



5 Gauss-Manin connection 

Let S{t) G C[t] such that 

n+l 

S{f)dx = dfA rjf, 7]f = Y,{-^^^Pidx^ G ^r-^ 

i=l 

For instance one can take S{t) := det{Af — t.I), where Af is the multiplication by / linear 
map form the Milnor vector space of / y := C[x]/iacob{f) to itself. This definition of S 
is implemented in the procedure S of brho.lib. The Gauss-Manin connection associated 
to the fibration / on H" turns out to be a map 

v.H"^ H'SMiPdx]) = i(Qp-Pf(f))d-\ p e CN 



where 

n+l 

^dxi^' ' dxi 



(11) «^ = E(|^«+^||) 

i=l 

satisfying the Leibniz rule, where for a set C C C by H'J^ we mean the localization of H" 

on the multiplicative subgroup of H" generated by t — c, c G C. Using the Leibniz rule 
one can extend V to a function frc 
times, makes sense. It is given by 



one can extend V to a function from H'^^ to itself and so the iteration V'^ = V o V • • • V k 



vyfc _ Vfc-1 ° Vfc-2 o • • • o Vc 
^^^^ ^ ~ 
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where 

Vfc : H" ^ H", Vk{[Pdx\) = [{Qp - {k + l)S'{t)P)dx] 
To calculate V : H' ^ H'q we use the fact that 

V'a; = — , uj(^H' 

df 

where d : H' ^ H" is taking differential and is well-defined. See §3 of ^ for more details 
on V. Usually the iteration of the Gauss-Manin connection produces polynomial forms 
with huge number of monomials. But fortunately our Brieskorn module H" (resp. H') has 
already the canonical basis ojp, £ I (resp. r/^, (3 € I) and after writing V the obtained 
coefficients are much more easier to read. The procedure nabla of brho . lib uses the 
formulas and (\12^ and computes V and its iterations. In H" one can write 

S{t)\/{uJf3) = ^ Pl3,l3'UJf3', Pf3,f3' G C[t], deg(p/3,/3') < deg(S') -1 + Ap-Af3' 
13' al 

The bound on degrees can be obtained as follows: 

Sujp = df Ar], ^ ds + dAp = d + deg(r/) 
This is because / is tame(see the proof of Lemma 4 of 

deg(p^,^/) < -Ap, = s-l + Ap-Afi, 

The procedure nablamat in brho. lib calculates the matrix g^^y The Gauss-Manin 
connection V has two nice properties: 

1. Griffiths transversality theorem: For alH = 1,2,..., n+\ we have S{t)'V[F^) C F^~^. 

2. Residue killer: For all lo £ H there exists a A; S N such that V^lu G Wn 

For the first one see 0. The second one for n = 1 is proved in Lemma 2.3 of "3]. The 
proof for n > 1 is similar and uses the fact that the residue as a function in t on a cycle 
around infinity is a polynomial in t. 



6 The numbers dp, /3 ^ I 

Let / be a tame polynomial with the last homogeneous part g, F be its homogenization 
and 

OF 

V = C[x, xo]/ < — I i = 1, 2, . . . , n + 1 > 

OXi 

We consider y as a C[xo]-module and it is shown in [Hj that V is freely generated by 
x^ := G /}. Let 

Ap-.V ^V, Af{G) = I^G, G £V 

OXq 

Proposition 3. The matrix of Ap in the basis x^ is of the form d. [j;q "''^'c^^^'], where 
Kjs^pi := d — 1 + deg(x^) — deg(x'^ ) and Af := is the multiplication by f in the 

Milnor vector space of f . In particular, if A pi — Ap > 1 then cp^pi = and 

det(AF - t.x^Q-^I) = dei{Af - t.I)x^Q'^^'' 
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Using the above Proposition, the procedure muldF calculates Ap. 

Proof. Since the polynomial F is weighted homogeneous, we have Yl7=o '^i^i§§^ — ^--^ 
dF 
dxQ 

n+1 

dXi 



and so = d.F in V (Note that ao = 1 by definition). Let 

(13) F.xf^ = ^ x'^'cp^p'ixo) + J2 '^/3,/3'(^o) € C[xo], Qi € C[xo,x] 



f3'Gl i=l 

Since the left hand side is homogeneous of degree d + deg(x^) we can assume that the 
pieces of the write hand side are also homogeneous of the same degree. This can be done 
by taking an arbitrary equation (|13|) and subtracting the unnecessary parts. □ 

Let C be a finite subset of C and C[t]^ be the localization of C[t] on its multiplicative 
subgroup generated hy t — c, c G C and Ft = F — I.Xq. From now on we work with 
C[t](j[xo,x] instead of C[2;o,a;] and redefine V using C[t](=i[xo, x]. Let 

dF, dF 

y^ = CM^[xo,x]/<— ,— ,|z = l,2,...,n + l> 

It is useful to reformulate V^, in the following way: Let R := C[t](^[3;o] be the set of 
polynomials in xq with coefficients in C[t]^ and At = Ap — t.d.XQ~^F We have 



V^ = V/ < —^q \qeV>= R^'/At.R^ 

OXo 



Here R^ is the set of /u x 1 matrices with entries in R. We consider the statement: 

*(C): There is a function /3€/^(i^€NU {0} such that the C[t] (^-module V^j is 
freely generated by 

(14) {x^°x^,0</?o <d^-l,/3G/} 

To prove the statement *(C) we may introduce a kind of Gaussian elimination in At and 
simplify it. For this reason we introduce the operation GE{Pi, (32, Ps). For /? G / let {At)/3 
be the /3-th row of At. 

• Input: At, Pi, (32, (33 G / with Ap^ < Ap^. 

Output: a matrix A[ and a finite subset i? of C. 
We replace {At)f3^ with 

(^t)/32,/33 



(A 



t)l3i,f32 



and we set B = zero{c{t)), where (^t)/3i,/32 = c(t).x^''^''^^ . Since for all /34 G / we 
have 

The obtained matrix A[ is of the form [xq '''''' ^/] and c'^^ = 0. If the matrix Bt is 
obtained from At by applying the above operation and B C C then At.R^ = BtR^. 

We give an example of algorithm which calculates d^'s for for some finite set C C C: 
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• Input: At 

Output: dp, (3 & I and a finite set (7 C C 

We identify / with {1, 2, . . . , //} and assume that 

Pl<P2^ Af3, > Ap^ 

The algorithm has /U steps indexed by /? = /x, — 1, . . . , 1. We define the set C to 
be empty. In /3 = ^ we have A{(3) = At. In the step j3 we find the first /3i such 
that A{(3)p^p^ / and put dp^=d-l + deg(x^) - deg(x^i). For /?2 = /? - 1, • • • , 1 
we make GE{(3, P2, Pi) and define C = C U U^^^-^^Bp^, where Bp^ is obtained during 
GE(P, P2, Pi)- The numbers dp^s obtained in this way proves the statement *(C). 

The advantage of this algorithm is that in many cases it gives C = C. We do not have a 
proof for *(C). One can also fix a value c G C\C and apply the above algorithm for Ac. 
In this case we do note care about C during the algorithm. The obtained dp^s make the 
statement *{C) true for some C C C with c ^ C. This algorithm is implemented in the 
procedure dbeta of brho.lib. We prove the following weak statements: 

Proposition 4. There is a function P G I ^ dp G NU {0} such that the C[t]c-niodule V' 
is generated by {xq°x^,0 < Po < dp — 1, P & I}. 

Proof. We have 

V = R>'/AtR^' ^ At^R^'/R^' = \, ,J R^ 



The isomorphism b in the middle is obtained by acting At from left on R'^ and adj makes 
the adjoint of a matrix. Now for P £ I let dp be the pole order of /3-th arrow of ^(d-i) ■ 

The numbers dp are the desired numbers. It is easy to see that {x^^xf^, < Po <dp,P e 1} 
generates V . □ 

Proposition 5. There is a subset C C C such that the statement *{C) is true with 
dp = d - 1,P e L 

Proof. We identify / with {1, 2, . . . , ^} and assume that 

Pi<p2^ Ap, > Ap^ 

By various use of operation GE on At we make all the entries of {At)p^^ = 0, /3 G 
We repeat this for {At)p^^-i = 0,/3 G — 1} and after //-times we get a lower 

triangular matrix. We always divide on a polynomial on t with leading coefficient one and 
so division by zero does not occur. □ 

Proposition 6. Let *{C) is valid with dp,P G /. Then 

Ap <n + l, dp < d{n + 2 - Ap), ^dp = fi{d - 1) 

Proof. The first one is already in Steenbrink's Theorem ^ The second inequality is ob- 
tained by applying the first inequality associated to F — cxq for some c G C\C: 

dp -1 + 1 

= ^^ + — < ^ + 2 

The Milnor number of F — cxq is X^^g/ dp and equals to the Milnor number ol g — cxq 
which is — 1). □ 
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7 Main theorem of fS] 

Suppose that *{C) is valid with dp, (3 G I. Define 



rk 

-'n+l 



{(3 e I \ Ap = n + 1 - k}, = {(3 e I \ Af3 + ^ < n + 1 - k < + ^} 



Theorem 2. For a tame polynomial f , the set 
with 



B'n+i = {V"-'^^/5 I P G i?^ = {V--'u;p I /? € I^}, 

is a basis of H'^ compatible with the mixed Hodge structure. The same is true for H'^ 
replacing V^~^ujji with V'^'^^'^rjp. 

Unfortunately, this theorem gives us a basis of a localization H compatible with mixed 
Hodge structure. In iilUI we have computed such bases for the Brieskorn module itself. 

To handle easier the pieces of the mixed Hodge structure of we make the following 
table. 







1 



jn jn 

-'n -'n+l 



2 

T-n — L 



n 



n + l 



rn— 1 jn—l Tn—2 
-'n+l -'n 



n 

-'n+l 



n 



The procedure Imk of brho . lib gives us x^, /3 G m 



order /. 



n jn—l 

w n ' ' 



rn— 1 



' -'n' -'n+l' -'n+l ) • • • ) -'n+l' 



^, m = n,n + 1, k = 0,1, ... n with the 
In the case n = 1 we have the table 



1^ — W 

-'1 -'2 



ll = {pel\Ap + ^<i<Ap + ^}, l',={PeI\Ap + ^<2<Ap + ^} 

4 = {(3eI\Ap = 1} 

The forms ujfs, P ^ ll form a basis of n Wi and the forms W/j,/? € form a basis of 
H"/Wi. Now to obtain a basis of Tyi/(-F^ n Wi) we must modify Vuj/s, (5 G 1°. 

The procedure changebase calculates the matrix of the basis of the Brieskorn module 
H'~, obtained in Theorem[7|in the canonical basis ijJf},j3 € I. 



8 Picard-Fuchs equations 

It is a well-known fact that for a polynomial / € C[x] and lo ^ H the integral I{t) := /^^ 
satisfies 

(15) CEP^it)^)It = ^, Pi{t)^C[t] 

called Picard-Fuchs equation, where 5t £ Hn{Lt,1,) is a continuous family of topological 
cycles. When / is tame, it is possible to calculate pi as follows: 
We write 

^\^) = '^Pi,fS^P 
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and define the k x fj, matrix A = [pj^/j], where i runs through 1,2, ... ,k and [3^1. Let 
k be the smahest number such that the the rows of ^fc_i are C(i)-hnear independent. 
Now, the rows of Ap^ are C(t)-hnear dependent and this gives us (after multiphcation by 
a suitable element of C[t]) 

k 

^p,(t)V^(a;) = 0, Pi{t)eC[t] 

i=0 

Using the formula and integrating the above equality, we get the equation (|T^ . The 
procedure PFeq from the library brho.lib calculates pj's in ()15|1. 



9 Polynomials in one variable, n = 

The theory developed in f^ldoes not work for the case n = 0. For a polynomial of degree 
d in one variable dim.(H^{Lt,C)) = d but fi = d — 1. However, if we use the following 
definition of homology and cohomology for a discrete topological space M, 

Hq{M, Z) = {m = ainii \ai ^TL, rrii & M \ deg{m) = Oj = 0} 

i i 

H^{M,C) = {/ : Ho{M,Z) C linear}/{/ [ f is constant on M} 

then 

H' = C[x]/C[f], H" = C[x]dx/f'C[f]dx, I = {l,x,x'^ , . . . ,x'^-^}, fi = d - I 
In this case 

uj = ^aiuj{pi), where 6 = ^aiPi, ai £ Z, pi e f^^{t), uj £ H' 
-'^ i i 

If, for instance, /' = has d distinct root then every vanishing cycle in Lt is a difference 
of two points of Lt. The set B = {x, x^ , . . . , x'^~^^ form a basis of W and its V which is 
{dx, xdx, . . . x'^~'^dx\ (up to multiplication by some constants) form a basis of B." . The first 
fact is easy to see. We write / = asx!^ + /o and for a polynomial p{x) € C[x] whenever we 
find some x'^ we replace it with ^~J^ and at the end we get p{x) = po{f) + Yli=i Piif)^^ or 

equivalently p = '}2i=i Pii^)^^ ^' ■ There is no C[t]-linear relation between the elements 
of B because B restricted to each regular fiber is of dimension d. We write 



d-2 d-2 



p{x)dx = Tq,{f)x'dx + qa-i{f)x''-'dx = (V%(/)x^dx - ^^^dx) + '^-^^^^dx 

and this proves the statement for H" . 

The proposition Q) can be stated in the case n = as follows: The only case in which 
dAp — = is when n = and P = 1. In the case n = for P 7^ 1 we have 

^■Pdx = — ^gdx 

OXi d.Ap — Ui oxi 
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and if P = 1 then -^.Pdx is zero in H" . The argument in Q and (jSJ can be done also 
in the case n = 0. In this case if 

(16) ri = ^Pf5{g)r}p+p{g), p,Pf3 G C[t] 

where each piece in the right hand side of the above equahty has degree less than deg(ry) 
then 

(17) dri = ^{pp{g)Ap + pp{g)g)up +p'{g)dg 

Based on this observation, the procedure linear, linearp, works for the case n = 0. 

In the case n = 0, we have only the set Jg = {Ajs + ^ < 1 < Ajs + -^} and this is equal 
to /. We have dp < d.{n + 2 - Ajs) = 2d- [5 -1 = and Ap = We conclude that 

d<dp + (3 + l<2d 

Now the infinitesimal Hilbert problem (see ^ Problem 7) can be stated in the case n = 0. 
Can one give an effective solution to this problem in this case? The positive answer to 
this question may give light into the the problem in the case n = 1. 

10 Examples 



For all the examples bellow we run 

>LIB "brho.lib"; 
>LIB "matrix, lib" ; 



10.1 Examples, n = 

For examples of this section we run 
>ring rO=(0,t) ,x, dp; 

Example 1. / = a;'^ - 5a:, P = {e* | i = 0, 1, 2, 3}, C = {-4e' 
i — 0, 1, 2, 3}, where e — e ^~ 1 is the d-th root of unity. 

> int d=5; poly f=x"d-d*x; okbase CstdC jacobCf ) ) ) ; 
.[i:-x3 

.[2>x2 

.[3]-x 

.[4]-l 

> AbetaCf ) ; 
.[1,1] -4/5 
.[2,1] =3/5 
.[3,1] -2/5 
.[4,1] -1/5 

> poly Sf-S(f); Sf; 
(14-256) 

> list l=nablamatCf ,Sf ) ; 

> 1[1] ; print (1 [2]); 
l/(5t4-1280) 

(-t3), 128, (-48t) , (16t2) , 
C4t2), (-213), 192, (-64t), 
(-16t) , (8t2) , (-3t3),256, 
64, (-32t) , (12t2) , (-4t3) 

//This is the matrix of nabla in the canonical 
//basis x"3,x"2,x"l,l. 

>PFeq(f ,1) ; 
.[1,1] -6144 
.[l,2]=(35625t) 
.[l,3]-(33376t2) 
.[l,4]=(8750t3) 
.[1,6]- (626t4- 160000) 



The residues of -r3t poles satisfy the Pieard-Fuehs equation 



6144 + 356254 h 33375t^ 8750t^ h 

dt 9*2 St3 



(6254* - 160000) — - = 

10.2 Examples n = 1 

For the examples bellow we define 

ring rl=(0,t), (x,y) , dp; 
Example 2. / — xy(x y — 1). 

> poly f- x2y+xy2-xy ; 

> poly g-lasthomo(f ) ; g; 
x2y+xy2 

> okbase (std(jacob(g) )) ; 
-[l]=y2 

-[2]-y 
_[3].x 
_[4].l 

> print(muldF(f-par(l))) ; 
(-3t+l/18)»x2,-l/18*x3, 0, 0, 
l/6*x, (-3t-l/6)»x2,0, 0, 
l/6*x, -l/6*x2, (-3t)*x2,0, 

1/2, -l/2»x, 0, (-3t)»x2 

> poly Sf-S(f ) ; Sf ; 
(t4+l/27t3) 

//We can take Sf =t* (t+1/27) ; 

> list ll=nablainat(f ,Sf) : 

> 11 [1] ; " "; print (11 [2]); 
l/(54t2+2t) 

(18t+l) , (-18t-l) ,0, (-2t) , 
1, -1, 0,(-6t), 
1, -1, 0,(-6t), 
3, -3, 0,(-18t) 
// 

> dbeta(f ,par(l)) ; 
0,2,2,4 

> Imk(f ,par(l)) ; 
[«: 
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[1]: 
[1]: 

1 

[2]: 
[1]: 

1 

[2]: 
[1]: 
[1]: 

X 

[2]: 

y 

> list 13=changebaseCf ,Sf ,paxCl)) ; 

> print (13 [1]); " print (13 [2]); det(13[2]); 
l,3/(54t2+2t),l,l 

0. 0, 0,1, 

1, -1,0, (-6t), 
0,0, 1,0, 
0,1, 0,0 
1 

// 

> dbeta(f ) ; 
2,2,2,2 

> Iiiik(f); 
[1]: 

[1]: 
[1]: 
1 

[2]: 
[1]: 

y2 

[2]: 
[1]: 
[1]: 

X 

[2]: 

y 

> list 12=changebase(f ,Sf ) ; 

> print (12 [1] ) ; " "; print (12 [2] ) ; det(12[2]); 
l,l/(64t2+2t),l,l 

0, 0, 0,1, 

(18t+l) , (-18t-l) ,0, (-2t) , 
0, 0, 1,0, 

0, 1, 0,0 

(18t+l) 

//The obtained basis does not work for the fiber c=-l/18 
// 

> PFeq(f,l, Sf); 
.[1,1] =6 
.[l,2]=(54t+l) 
.[l,3]=(27t2+t) 
.[1,4]=0 
.[1,5]=0 

We get the following basis of h" compatible with mixed Hodge 
structure. 



/ = xy{x + y — \) 





[1] 


Gr" Grf^H" 


[y'] - [y] - 6t[i] 


Gr\,Gr^H" 


N,[y] 



The integrals / = /^^ 'T—^ satisfy the Picard-Fuchs equation 

dl o d^I 

6 + (54t + 1)— + (27*'' + t)— = 
dt dt-^ 

Example 3. / = 2{x^ + y^) - 3(x^ + y^), P = 
{{0, 0), (0, 1), (1, 0), (1, 1)}, C = {0, -1, -1, -2}, 

> poly f= 2«3+2*y3-3*x2-3*y2 ; 

> poly g=lasthomo(f ) ; g; 

2*x3+2*y3 

> okbase (std( jacob(g) ) ) ; 
.[l]-xy 

.[2]=y 
.[3]=jt 
^M-l 

>S(f); 
(t4+4t3+5t2+2t) 
//We can put 

>poly Sf-t*(t+l)*(t+2) ; 

> list 12=changebase(f ,Sf) : 

> print (12 [1] ) ; " "; print (12 [2] ) ; det(12[2]); 
l,-l/(6t+12),l,l 

0, 0,0,1, 
-2,1,1,0, 
0, 0,1,0, 
0, 1,0,0 
-2 



pi.,/ J n 


[11 


Gr^pGrf if" 


I2xy — X — y] 


Gr^Gr^ H" 





Example 4. / — + y^ — x. 



> poly f= x4+y4-x ; 

> poly g=lasthomo(f ) ; 

> okbase(stdCjacob(g) ) ) ; 
_[l]=x2y2 

-[2]=xy2 

.[3]=x2y 

.[4]=y2 

.[5]=xy 

_[6>x2 

.[7]-y 

.[8].x 

.[9].l 

> poly Sf-S(f) ; Sf ; 
(t9+81/256t6+21S7/65536t3+19683/16777216) 
//We can take 

>Sf-t-3+27/256; 

> dbetaCf ,parCl)) ; 
2,2,2,6,2,2,6,2,6 

> Iiiik(f ,par(l))i 
[1]: 

[1]: 
[1]: 
1 

[2]: 

X 

[3]: 

y 

[2]: 
[1]: 

y 

[2]: 

y2 

[3]: 
x2y2 

[2]: 
[1]: 
[1]: 

x2 
[2] : 

xy 
[3]: 

y2 

> list 13=changebase(f ,Sf ,par(l)) ; 

> print(13[l]) ; " "; print (13 [2] ) ; det(13[2]); 
1,1, 1 ,4/(256t3+27) ,24/ (256t3+27) , 1/ (256t3+27) ,1,1,1 



0, 


0, 


0,0,0,0,0, 0,1 


0, 


0, 


0,0,0,0,0, 1,0 


0, 


0, 


0,0,0,0,1, 0,0 


0, 


0, 


9,0,0,0,(-16t2),0,0 


3, 


(-2t) 


,0,0,0,0,0, 0,0 


(128t2) 


,9, 


0,0,0,0,0, 0,0 


0. 


0, 


0,0,0,1,0, 0,0 


0, 


0, 


0,0,1,0,0, 0,0 


0, 


0, 


0,1,0,0,0, 0,0 



(2304t3+243) // 9»256«Sf; 

> matrix A-13E2J ; 

> A [6,1. .ncolsCA)]- 

(C-128*t2)/3)*submat CA,5,1 . .ncols CA) )+submat (A ,6 , 1 . .ncols(A)) ; 

> A[5,l. .ncols(A)]- 

2*t*submat (A, 6,1. .ncols(A) )+submat (A,5, 1 . . ncols (A) ) ; print (A) ; 
0,0,0,0,0,0,0, 0,1, 
0,0,0,0,0,0,0, 1,0, 
0,0,0,0,0,0,1, 0,0, 
0,0,9,0,0,0, (-16t2) ,0,0, 
1,0,0,0,0,0,0, 0,0, 
0,1,0,0,0,0,0, 0,0, 
0,0,0,0,0,1,0, 0,0, 
0,0,0,0,1,0,0, 0,0, 
0,0,0,1,0,0,0, 0,0 

We obtain the following table 



r 4.4 

J = X + y — X 



Gr|,Grf H" 


[1], W,M 


GrlGr't h" 


9[x-'y] - let^fo], \x^y^],[xy^] 


GrJ^Gr^' H" 


[^'%[^y\.[y^] 



We make the following remark 
> reduce(9*x2*y-16*Cf "2)*y, stdCjacobCf ))) ; 
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